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RAPID  GEOACOUSTIC  INVERSION  METHODS  USED 
IN  LWAD98-2  SEA  TEST 


1)  Introduction 

The  characterization  of  bottom  geoacoustic  and  statistical  properties  is  required  for 
accurate  prediction  of  sonar  performance  in  littoral  waters.  In  many  of  the  areas 
considered  to  be  sites  of  future  Navy  operations  little  is  known  of  the  bottom 
composition,  particularly  at  depths  more  than  a  meter  into  the  bottom.  The  real-time 
capability  of  remotely  estimating  seabed  geoacoustic  and  statistical  properties  is 
especially  important  for  Navy  vessels  operating  in  forward  areas  over  previously 
unsurveyed  bottoms. 

During  the  Littoral  Warfare  Advanced  Development  (LWAD)  98-2  sea  test,  the  concept 
of  rapid  environmental  characterization  using  broad-band  pulses  was  tested  off  the 
southwest  coast  of  Florida  in  February/March  1998.  The  experiment  area  was  located  in 
the  Gulf  of  Mexico  near  the  Southwest  Florida  shelf  break  (see  Figure  1).  In  addition  to 
NRL  single-frequency  (1.5,  2.0,  2.5,  3.0,  3.5,  4.0,  and  4.5  kHz)  monostatic  scattering 
measurements,  broad-band  (1.5  to  4.5  kHz)  pulses  were  used  to  invert  geoacoustic  and 
statistical  properties  at  site  SR2  (25  36.26N,  84  07.95W).  The  idea  of  retransmitting  the 
portions  of  the  scattered  signal  for  the  enhancement  of  signal  returned  from 
corresponding  annular  patches  is  also  tested.  This  report  briefly  describes  the  inversion 
techniques  and  their  usage  in  a  series  of  proof-of-concept  and  engineering  trials  during 
the  LWAD98-2  experiment.  The  next  two  sections  describe  a  bottom  backscattering 
model  and  a  normal-incidence  reflection  model  which  were  used  for  inverting 
geoacoustic  and  statistical  properties  from  backscatter  data.  Experiment,  data  analysis, 
and  results  are  discussed  in  Section  4.  Finally  a  brief  summary  is  given  in  Section  5. 

2)  Jackson/Turgut  bottom  backscattering  model 

Bottom  backscattering  strength,  (0)  is  defined  as  the  decibel  equivalent  of  the  sum  of 
interface  and  volume  scattering  cross  sections  (Urick,  1983) 

S,(0)  =  lOlogiob,(0)  +  <T,(0)]  (1) 

where  e  is  the  grazing  angle,  (0)  and  cr^  (0 )  are  dimensionless  backscattering  cross 
sections  per  unit  solid  angle  per  unit  area  due  to  interface  roughness  and  volume 
scattering,  respectively.  Both  cross  sections  are  calculated  using  perturbation 
approximations,  and  corrections  for  shadowing  and  large-scale  bottom  slope  are  included 
in  the  context  of  the  composite  roughness  approximation  (McDaniel  and  Gorman,  1983) 
{6)  =  S  (0,  s)  Fid,  ,  5) 

(r,(0)  =  5(0,5)F(0,(j^,,5) 

where  5(0)  is  the  shadowing  correction,  F(6,g,s)  is  a  slope  averaging  integral,  and  s  is 
the  rms  slope.  A  detailed  evaluation  of  shadowing  correction  and  rms  slope  averaging 
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integral  can  be  found  in  the  APL-UW  Technical  Report-9407  (1994).  The  bottom 
backscattering  cross-section  Op^iO)  in  the  Rayleigh-Rice  approximation  is  given  as 

(Kuo,  1963) 

(T^^(0)  =4k'^  sin^  e\Yiefw{2k' cosd)  (3) 


where  k’  is  the  acoustic  wavenumber  in  the  water.  The  complex  function 7(0)  is  defined 


Y{G)  = 


(^-l)^COS0  +  ^^  -t' 


qsmd  +  {'i 


-cos^  0)* 


where  q  =  Pq  /  is  the  ratio  of  the  density  of  the  sediment  to  the  density  of  the  water, 

t  =  (1  +  /5)/v  is  the  ratio  of  the  complex  wavenumber  in  the  sediment  to  the  real 
wavenumber  in  the  water,  8  is  the  loss  parameter,  and  v  =  Cq  I  Cf  is  the  ratio  of  the 

sediment  sound  velocity  to  water  sound  velocity.  In  the  volume  scattering  model  of 
Turgut  (1998),  after  including  refraction  and  transmission  loss  at  the  interface,  volume 
scattering  strength  as  a  function  of  grazing  angle  is  written  as 


o-,v(e)  = 


2sin0, 


where  r\  =  [4;rY  +  Im(fe)]~'  is  the  inverse  attenuation  coefficient,  6^  is  the  refracted  angle 
in  the  sediment,  and  Y  is  the  volume  scattering  cross  section.  Two-way  transmissivity 
r(0)  in  the  above  equation  is  defined  as 

4^vsm»sin9,  ^  (6) 

(^vsin0  +  sin0,.) 

Considering  a  transversely  isotropic  structure  of  sediment  with  sound-speed 
inhomogeneities  having  von  Karman  type  of  spectrum,  the  resulting  volume  scattering 
cross  section  can  be  written  as 

k  a  bV  —  r(m) — cos^  0^ -f-0^  sin^  0j]  ,  m>3/2,  (7) 

24k  r(m-3/2)L  V  n 

where  a  where  is  the  horizontal  scale  factor,  b  is  the  vertical  scale  factor,  T  is  the 
Gamma  function,  and  m  is  the  spectral  exponent.  The  parameter  is  defined  as 

V  (8) 

^  y  4.y  _2 

where  Yp  =  Pr  /  Po  ratio  of  the  grain  density  to  the  bulk  sediment  density  and 

Y^  I  kTq  is  the  ratio  of  the  grain  compressibility  to  the  bulk  sediment  compressibility. 


3)  Inversion  of  seabed  geoacoustic  properties  form  reflection  data 

As  a  forward  model,  the  seabed  is  modeled  as  porous  homogenous  layers  and  the 
wavefield  is  conveniently  decomposed  into  upgoing  and  downgoing  waves  within  each 
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homogeneous  layer.  A  matrix  equation  relates  the  upgoing  and  downgoing  waves  at  the 
bottom  interface  of  layer  i  to  the  upgoing  and  downgoing  waves  at  the  top  of  layer  /+!  as 


(  Ri  ^ 

T<  —!_J.  _L 

fu  ^ 

^i+I 

'  T  T 

1  1 

R\  1 

T,  TJ 

A+l  ) 

where  Ri,T{  are  the  reflection  and  transmission  coefficients  at  interface  i  for  waves  that  are 
incident  above  and  7”i  are  the  reflection  and  transmission  coefficients  for  waves  that 
are  incident  from  below.  The  details  of  the  derivation  of  seiosmogram  equations  for  a 
layered  poro- viscoelastic  medium  can  be  found  in  Turgut  and  Yamamoto  (1988). 

A  previously  developed  time-domain  inversion  algorithm  (Turgut  and  Wolf,  1996)  was 
used  for  inverting  seabed  geoacoustic  properties  from  normal-incidence  reflection  data. 
First,  envelope  functions  of  the  matched-filtered  acoustic  returns  are  calculated  by  taking 
the  absolute  values  of  the  acoustic  returns  in  the  form  of  analytic  signals.  The  matched 
filtering  achieves  a  signal  processing  gain  over  the  background  noise  and  provides 
compressed  matched  filter  output.  Then,  instantaneous  frequencies  are  calculated  for  a 
crude  estimation  of  attenuation  profile  by  comparing  the  center  frequency  shift  in  the 
acoustie  return  with  that  predicted  by  using  the  poro-viscoelastic  attenuation  model. 
Assuming  an  average  sound  speed  value  (1500  m/s)  in  the  water  column,  water  depth  is 
calculated  from  the  delay  time  between  the  transmission  peak  and  the  first  return  peak. 
Next,  porosity,  shear  modulus,  and  permeability  are  inverted  for  the  first  sediment  layer 
by  using  a  global  search  algorithm  within  a  model  space.  The  model  space  is  further 
constrained  by  well-established  semi-empirical  relations  among  the  model  parameters 
such  as  porosity,  shear  modulus,  and  attenuation.  Random  perturbations  are  applied  to 
corresponding  model  parameters  to  minimize  an  objective  function  (quadratic  deviation 
between  the  data  and  calculated  synthetic  seismogram  for  the  first  reflector).  Acoustic 
parameters  such  as  compressional  speed  and  density  are  also  ealculated  from  the  inverted 
model  parameters  using  the  poro-viscoelastic  theory. 

4)  Experiment  and  Data  Analysis 

a)  Inversion  Using  backscatter  data 

The  monostatic  scattering  data  were  collected  by  using  an  existing  NRL  bottom- 
backscattering  system.  The  NRL  bottom  backscattering  system  uses  a  USRD-G81 
transducer  as  a  source  (source  level  187  dB  re  1  |xPa  @  Im)  and  two  16-element  hardened 
hydrophone  arrays  with  0.15  m  element  spacing.  Prior  to  the  experiment,  several 
wideband  (1.5  kHz  to  4.5  kHz)  50  ms  duration  LFM  waveforms  were  designed  and 
digitally  stored  on  a  personal  computer  equipped  with  a  digital-to-analog  conversion 
board.  A  transmission  schedule  was  designed  for  each  site  based  on  the  source/receiver 
depth  and  water  depth.  Typically,  50  ms  10-25%  cosine  tapered  waveforms  were 
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transmitted  every  15  s,  allowing  a  sufficient  decay  of  the  reverberation  tail.  The 
backscattering  strength  was  calculated  by  using  the  relation  =  101og,o  where  the 
scattering  coefficient  is  given  as  (Bunchuk  and  Zhitkovskii,  1980) 


M,= 


R 


lo 


(10) 


V  VY 

In  the  above  equation,  R  is  the  distance  from  source  to  the  insonified  area  at  the  water- 
sediment  interface,  is  the  sound  speed  in  water,  /q  is  the  intensity  of  the  transmitted 


signal  at  a  distance  Rq  from  the  source,  T  is  the  finite  duration  of  the  transmitted  pulse, 
and  is  the  intensity  of  the  scattered  signal.  Since  the  transmitted  and  scattered  signals 
were  recorded  by  the  same  hydrophone  array  without  clipping,  relative  levels  of  7^  and 
/q  were  used  in  the  backscattering  calculations.  In  Figure  2,  the  backscatter  levels  are 
depicted  as  a  function  of  time  and  steering  angle,  calculated  from  the  beam-time  data  by 
using  the  above  equation.  The  bottom  and  surface  scattering  begin  at  about  0.1  s  with 
bottom  scattering  component  being  dominant  due  to  a  null  in  the  upward  direction  of  the 
source  beam  pattern.  The  backscatter  levels  for  frequencies  with  0.5  kHz  windows  are 
also  shown  in  Figure  2.  The  comparable  bottom  scattering  levels  for  each  frequency  band 
indicate  minor  frequency  dependency  for  this  site.  The  backscattering  levels  were  also 
comparable  with  those  of  the  recently  reported  single-frequency  backscattering  scattering 
analysis  (Kunz,  1998).  Table  1  shows  the  geoacoustical  and  statistical  parameters  of  the 
site  SR2  area  estimated  from  the  backscattering  data  at  3  kHz.  The  corresponding  a 
posteriori  probability  density  (PPD)  functions  are  plotted  in  Figure  3  with  the  mean  and 
one  standard  deviation  of  the  estimation.  PPDs  are  good  indicators  of  importance  and 
uniqueness  of  each  parameter.  The  estimated  geaoacoustic  and  statistical  parameters 
indicate  the  presence  of  t}q)ical  sandy  bottom  in  the  area.  Figure  4  shows  the  measured 
backscattering  data  and  the  model  prediction  using  the  inverted  parameters.  Notice  that 
the  model  predictions  are  usually  within  an  assumed  3-dB  measurement  error. 


b)  Inversion  of  normal-incidence  reflection  data 

The  seabed  geoacoustic  properties  such  as  compressional  speed,  density,  and  attenuation 
are  also  estimated  by  using  matched-filtered  and  endfire-beamformed  data.  Figure  5 
shows  matched  filter  output  of  direct  and  bottom  reflected  arrivals  by  steering  the  array  in 
both  end-fire  directions.  Notice  that  amplitude  of  direct  signal  is  about  four  times  larger 
than  that  of  first  bottom  reflected  signal  (after  the  removal  of  geometrical  spreading).  A 
strong  second  signal  right  after  the  bottom  return  indicates  a  2-3  m  thick  sediment  layer 
overlying  a  hard  bottom  (strong  reflector).  Figure  6  shows  the  instantaneous  frequency 
and  envelope  function  of  the  bottom  reflected  signal.  The  average  center-frequency  shift 
calculated  from  the  instantaneous  frequency  is  about  13  Hz/m  which  corresponds  to 
sandy  bottom  in  the  Biot/Turgut  attenuation  model  (Turgut,  1998).  For  comparison, 
model  prediction  of  the  center  frequency  shift  is  given  in  Table  2  for  different  types  of 
sediments.  Figure  7  shows  inverted  compressional  speed  and  density  profiles  for  the  site 
SR2  location.  The  values  are  comparable  to  those  of  a  geoacoustic  model  reported  in 
Bucca  and  Fulford  (1997).  However,  additional  “ground  truth”  data  such  as  sediment  core 
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measurements  are  needed  to  validate  the  inversion  method  applied  to  NRL  scattering 
system  data. 

c)  Time  reversal  in  bottom  scattering  measurements 

Application  of  time  reversal  in  the  bottom  backscattering  measurements  is  expected  to 
produce  robust  backscattering  strength  measurements  since  the  eorrelations  of  the 
retransmitted  signals  are  done  automatically  by  the  seabed  itself  (Clouet  and  Fouque, 
1997).  Thus,  a  reduction  in  the  measurement  errors  and  data  collection  at  lower  grazing 
angles  are  expected.  The  engineering  trial  of  the  time-reversal  software/hardware  system 
was  satisfactory  except  for  the  insufficient  filtering  of  the  returned  signal.  In  the 
retransmitted  signals,  the  low-frequency  noise  levels  were  much  higher  then  the  bottom 
scattering  signals  preventing  a  comprehensive  testing  of  the  method  itself.  This 
deficiency  has  been  removed  by  using  better  band-pass  filters  and  the  system  (and  the 
method)  is  planned  to  be  tested  in  future  sea  trials. 


5)  Summary 

Rapid  environmental  characterization  of  the  sea  bottom/subbottom  using  the  NRL 
scattering  system  has  been  performed  during  the  LWAD98-2  seatest.  The  statistical 
parameters  of  a  Jackson/Turgut  bottom/subbottom  scattering  model  were  estimated  from 
the  backscatter  data.  Seabed  geoacoustic  properties  were  also  inverted  from  normal- 
incidence  acoustic  data.  Results  are  in  agreement  with  the  previously  reported 
geoacoustic  data  for  the  area  indicating  a  2-3  m  thick  sandy  layer  over  a  hard  subbottom 
(coral  rock).  The  engineering  trial  of  a  time-reversal  backscattering  system  was 
satisfactory  except  for  the  insufficient  filtering  of  the  low-frequency  noise.  This 
deficiency  has  been  removed  and  will  be  tested  in  a  future  sea  trial.  Also,  the  inversion 
methods  should  be  further  validated  by  comparing  the  results  with  high-resolution 
sediment  coring  data  from  an  area  where  the  bottom  composition  is  well  documented. 
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Surface  spectral  exponent 

(y) 

3.6245  (+/-0.2101) 

Surface  spectral  strength 

(P) 

0.0066  (+/-0.0010) 

Loss  parameter 

(5) 

0.0159  (+/-0.0016) 

Aspect  ratio 

(a) 

5.7949  (+/-0.3895) 

Vertical  scale  factor 

(b) 

0.2269  (+/-0.1322) 

Sound  speed  variance 

(ft) 

0.0170  (+/-0.0029) 

Sound  speed  ratio 

(V) 

1.1301  (+/-0.0115) 

Volume  spectral  exponent 

(m) 

1.7342  (+/-0.0997) 

Table  1.  Inverted  geoacoustic  and  statistical  parameters  from  3  kHz  backscatter  data.  One 
standard  deviation  of  the  estimation  are  given  in  parentheses. 


Frequency  shift  (Hz/m) 

S=0.02 

8  =  0.05 

5=0.1 

5  =  0.2 

Silty  Clay 

0.67 

0.88 

1.15 

1.69 

Clayey  Silt 

1.28 

1.48 

1.75 

2.43 

Silty  Sand 

8.54 

8.73 

9.10 

9.79 

Fine  Sand 

14.57 

14.81 

15.17 

15.95 

Table  2.  Center  Frequency  shift  for  a  1. 5-4.5  kHz  signal  estimated  by  Biot/Turgut 
sediment  attenuation  model  for  different  frame  loss  coefficient  (6). 
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Recent  Carbonate  Reef 
Unidentified 


Figure  1.  Location  and  bottom  composition  of  LWAD98  bottom  scattering  experiment 
site  SR2. 
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Figure  3.  A  posteriori  probability  density  functions  of  geoacoustic  and  statistical 
parameters  inverted  for  the  SR2  site.  The  mean  and  one  standard  deviation  of  the 
estimation  are  also  shown. 
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Figure  4.  Measured  backscattering  data  and  model  prediction  using  inverted  parameters. 
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Alton  Tursut 


Figure  5.  A  typical  matched  filter  output  of  endfire  beamformed  data 


Figure  6.  Instantaneous  frequency  (dotted  line)  and  envelope  function  (solid  line)  for  the 
average  attenuation  estimation.  Regression  line  indicates  a  ~13  Hz/m  center  frequency 
shift.  Dots  indicate  relatively  high  signal  strength  data  selected  for  regression  analysis. 
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Figure  7.  Estimated  compressional  wave  speed  and  density  profiles  at  the  SR2  site. 
Depth  unit  is  in  meters  below  seafloor. 


